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Abstract. We study the stability of unitary quantum dynamics of composite systems 
(for example: central system + environment) with respect to weak interaction between 
the two parts. Unified theoretical formalism is applied to study different physical 
situations: (i) coherence of a forward evolution as measured by purity of the reduced 
density matrix, (ii) stability of time evolution with respect to small coupling between 
subsystems, and (iii) Loschmidt echo measuring dynamical irreversibility. Stability has 
been measured either by fidelity of pure states of a composite system, or by the so-called 
reduced fidelity of reduced density matrices within a subsystem. Rigorous inequality 
among fidelity, reduced-fidelity and purity is proved and a linear response theory is 
developed expressing these three quantities in terms of time correlation functions of 
the generator of interaction. The qualitatively different cases of regular (integrable) or 
mixing (chaotic in the classical limit) dynamics in each of the subsystems are discussed 
in detail. Theoretical results are demonstrated and confirmed in a numerical example 
of two coupled kicked tops. 

PACS numbers: 03.65. Sq, 03.65.Yz, 05.45.Mt 
1. Introduction 

Recently we have witnessed a strong revival of interest in the theoretical questions 
related to decoherence, in particular due to an immense practical potential of the 
upcoming quantum information processing technology Jl||. In order to design quantum 
devices capable of coherent quantum manipulation, one has to be able to control and 
minimize the decoherence due to an unavoidable weak coupling between the system 
of the device and the environment. Traditionally, one uses idealized harmonic heat 
bath models of an environment, and very often also harmonic description of the central 
system (e.g. 0). However, the rate of decoherence, or quantum dissipation, may 
depend on the intrinsic dynamics of the central system, whether be chaotic or not. It 
has been argued by Zurek |3], Q that the rate of decoherence as characterized by the 
Von Neumann entropy of the reduced density matrix increases with the rate which 
is given by the classical phase space stretching rate (Lyapunov exponent), for the 
quantum state which is initially given as minimal uncertainty Gaussian wave-packet 
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and for sufficiently short times. Recently, Jalabert and Pastawski have studied the 
so-called quantum Loschmidt echo 0, or quantum fidelity, which may be treated as a 
measure of dynamical reversibility under a slight change in the Hamiltonian, and found 
a similar relation to classical stretching rate for short times as for decoherence. This 
reflects correspondence between classical and quantum evolution of wave packets up to 
the Ehrenfest time — log ft ||. Later on, numerous papers appeared addressing related 
issues H- 

On the other hand, we have shown |9|, |l(J that the behaviour of quantum fidelity 



may be completely different, if either the time is longer than the Ehrenfest time, or if the 
initial (pure) state is more complex, e.g. random. In general, the rate of fidelity decay 
is given by the time-integral of autocorrelation function of the perturbation operator, 
and this is bigger the less chaotic is the dynamics thus making the fidelity lower, and 
vice versa. 

In previous papers [§, [10], [11], [L^, |13[ we have studied fidelity and the so-called 



purity fidelity characterizing the stability of quantum evolution, or quantum echoes, 
to the perturbation of dynamics in a generally coupled composite system, i.e. where 
both perturbed and unperturbed systems were coupled. In a situation where we are 
interested only in the properties of a central subsystem, one usually does not have the 
above general situation but a more specific one. Namely, the coupling to the environment 
is usually unwanted and small and so in an ideal (unperturbed) case we would like to 
have two decoupled systems. In that case we are interested in how the coupling of a 
central system to the environment changes the properties of the central system. So 
rather than comparing evolution in two general systems, in which both have coupling 
between subsystems, we have a situation where in an unperturbed case the systems 
are uncoupled and become coupled only because of the perturbation. This situation is 
studied in the present paper. 

We develop a unified theoretical framework to deal with different physical situations: 
(i) coherence of forward evolution as measured by purity of the reduced density matrix 
traced over the subsystem, (ii) stability of time evolution with respect to small coupling 
between the subsystems, and (iii) Loschmidt echo measuring dynamical irreversibility. 
Stability is measured either by fidelity of pure states of a composite system, or by 
the so-called reduced fidelity of reduced density matrices within a subsystem. We find 
general linear response formula expressing the fidelity, the reduced fidelity and the purity, 
in terms of time-correlation functions of the generator of the perturbation within the 
subsystems. We emphasize that the decay rates as given by linear response formalism 
are usually valid also in the regime of small fidelity /purity. Our general qualitative 
conclusion is that all the three quantities decrease slower with the increasing chaoticity 
of the dynamics in the subsystems. Our theoretical results are clearly demonstrated on 



a system of two coupled quantized kicked tops f!M . In addition we find some intriguing 



numerical results on algebraic long-time tails of some of the stability measures. 
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2. Characterizing stability and coherence of reduced time evolution 

Our system consists of a central system and an environment, henceforth denoted by 
subscripts "s" and "e", respectively The Hilbert space of a composite system is a 
direct product H = H s <8> TL e - We will compare quantum evolutions generated by two 
Hamiltonians, the unperturbed H, and the perturbed H$, 

H = H s ® l c + l s ® H e , H 5 = H + 5-V, (1) 

where H Sye acts only on the corresponding subspace 7i S)e and 5 is a dimensionless coupling 
strength, and V is a general perturbation operator which couples both systems. The 
usual measure of stability of overall unitary evolution on the total Hilbert space TC is 
the pure state fidelity (equivalent to quantum Loschmidt echo), which is the overlap 
between perturbed \4>s(t)) = ^(^^(O)) and unperturbed \i[>(t)) = U(t)\ip(0)) time 
evolving states, where U(t) = exp(—iHt/h) and Us(t) = exp(—iHst/h) are the unitary 
propagators. It turns useful to define an echo operator M$(t) [[U], |] 

M s (t) :=U\t)U s (t). (2) 

Writing the density matrices, 



Pit) 
Ps(t) 

PMs{t) 



\m)(m\ = u(t)p(o)uHt), (3) 

\Mt))(Mt)\ = Us(t) P (0)Ut(t), (4) 
M 3 (t)p(0)M}(t), (5) 

the fidelity can be concisely written as 

F(t) := \(Mt)\m)\ 2 = tr[p(0) PMS (t)]. (6) 

The echo operator (|2|) can be rewritten as 

M s (t) = texp(-iZ{t)6/ti), with E(t) := /' V(r)dr, (7) 

Jo 

where T is a time-ordering operator, and V(t) := W(t)VU(t) is the perturbation in the 
interaction picture. This representation of the echo operator ([FD is very convenient as 
it can be used |12| to derive various results on the behaviour of fidelity (|^) and related 
functions. As Pm<5(^) is nothing else but the density operator of the total system in the 
interaction picture, it satisfies 
d \S 

j t Pus(t) = jr[pus(t),V(t)], (8) 

where [A, B] := AB-BA. 

The fidelity is a measure of the distance within the whole Hilbert space 7i. On 
the other hand in the spirit of our study of a central system, we would like to have a 
quantity that would measure distance only within the central system space 7i s . So 
we should not care if the environmental states are corrupted by the perturbation, 

| Note that here the order of perturbed and unperturbed propagators is interchanged with respect to 
refs. jO, |l^] in order to facilitate exact partial tracing, but is otherwise inessential. 
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but would just like the evolution on the subspace H s to be preserved as closely 
as possible. With this aim we define a new quantity that we call reduced fidelity, 
which is a fidelity between the reduced perturbed p s s(t) := tr e ps(t) and unperturbed 
p s (t) : = tr e p(t) = U a (t)p s (Q)U%(t) density matrices, which start from the same product 
(disentangled) initial state \ip(0)) = |^(0)) s (g> \ip(0)) e and we use obvious notation 
U 8 (t) := exp (—iH s t/h). The reduced fidelity F R therefore reads 

F R (t) := tr s [p s (t)p sS (t)] = tr s [p 8 (0) tr e {pm(t)}]. (9) 

The reduced fidelity measures the distance between the two reduced density matrices. 
It can be interpreted either as an inner product between two reduced density operators 
propagated by two nearby hamiltonians, or as an inner product between the initial and 
the final reduced density operator after the echo dynamics. 

On the other hand, if we are interested only in the separability (disentanglement) 
of the final echo density matrix pMs(t), the relevant quantity is purity-fidelity JET], IT3] 

Fp(t) :=tr s [tr c p M 5(t)] 2 . (10) 

Now we shall use the fact that H is separable in two subsystems ([!]), so that purity- 
fidelity is in this case equal to purity of the coupled forward time evolution 

I(t) := tr s [tr c p 5 {t)f = tr s [p sS (t)} 2 . (11) 

In order to see that purity-fidelity is equal to purity in this situation we bring separable 
propagator (U — U s ® U e ) out of the innermost trace and use cyclic property of a trace 
opeartion in the definition of Fp(t) (|T0|) 

F P (t) = tr s [tr c (U^Usp(O)UjU)} 2 = tr s [£# tr c {U s p(0)Uj}U s } 2 = I(t). (12) 

This can be understood as a consequence of the purity being constant during the 
evolution with the separable hamiltonian H. However in the general case, the purity- 
fidelity is a property of echo dynamics and is different from the purity of forward 
dynamics. Since in this paper we are interested in the former case ([!]), we will from now 
on use a symbol I(t) instead of F P (t). 

Summarizing, all the three quantities, namely F(t), F-&(t) and I(t), measure the 
stability of a composite system to perturbations. The fidelity F(t) measures the stability 
of a whole state, the reduced fidelity gives the stability on H s subspace and the purity 
measures separability of pg(t). One expects that the fidelity is the most restrictive 
quantity of the three - pit) and p$(t) must be similar for F(t) to be high. For F R (t) to 
be high, only the reduced density matrices p s (t) = tr c [pit)] and p B s(t) = tr e [/>$(£)] must 
be similar, and finally, purity I(t) flTTD is high if only ps(t) factorizes. In a previous 
paper [|1]]] we proved an inequality F 2 (t) < Fp(t). Along the same line one can prove 
the following general inequality 

F 2 (t) < F*(t) < I(t). (13) 

Proof: Write p(0) = p s <S> p c . Uhlmann's theorem (noncontractivity of fidelity) J0J 
states for any pure states pus(t) an d p(0) that 

F(t) = tr [p(0)pMs(t)] < ti s [p s ti c {p M5 (t)}] = F K (t). (14) 
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Then, squaring and applying Cauchy-Schwartz inequality |tr(At.B)| 2 < ti(AA^)tr(BB^) 
we get F 2 (t) < F^(t) < tr s [tr e p M(5 (t)] 2 = F P (t) = I(t) with equality being satisfied only 
in the trivial case of Mg(t) = 1, i.e. when 5 = 0. 



3. Linear response 

Next we proceed by expanding all the three important quantities in powers of the 
perturbation strength 5. Although this procedure is very elementary, it greatly helps in 
understanding the behaviour of various measures of stability. The lowest order is given 
by various two-point time correlation functions of the perturbation and the decay time 
scale can be inferred from the behaviour of time integrals of this correlations. What 
is more, the dependence of this time scale on 8 and h as well as on the dynamics of a 
system (fast correlation decay or absence of correlation decay) is explicit. We start by 
expanding the echo operator 

M,(t) = l-^S(t)-^TS 2 (t) + --- (15) 
The leading order expansion of fidelity (|6|) is then 

F(t) = C(t), C(t) : = <£ 2 (t)> - (S(t)) 2 (16) 

where (•) denotes an expectation in the product initial state \i[>(0)) = 1 1 , 1 ) , with the 
general notation for a complete basis of Hilbert space \i,u) := \i) B <8> |i/) e . If explicitly 
written out, the coefficient C(t) is just an integral of autocorrelation function and reads 

C(t) = f f{(V{£)V{Q) - (V(0)(V(0)}dtd(. (17) 
Jo Jo 

Similarly, for the reduced fidelity F&(t) we obtain 

F R (t)= 1-(T) {C(t)-D(t)} (18) 
D(t) := (E(t)(Ps ® l.)E(t)> - (S(t)) 2 = ^ 1(1, H w, i)i 2 , 

and for the purity lit) 

I(t) = 1-2(|) {C(t)-D(t)-E(t)} (19) 
E(t) := (Y,(t)(l s <S> p e )E(t)) - (S(t)) 2 = l|E(t)|l, 1)| 2 . 

So far we have not specified any particular form of the perturbation V yet. To 
facilitate calculations, we now assume the simplest and physically well justified product 
form of the interaction 

V:=V B ® V c . (20) 

This is a very natural choice in the studies of decoherence. Such is the usual case 
where one writes V = x s (g) F e meaning the coupling of position times force. Henceforth, 
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the operator V e will be referred to as "force" and V s as "position". Following this 
assumption, the three coefficients C(t), D{t) and E(t) can be written out explicitly in 



terms of separate correlation functions over different spaces, namely («) f 



C(t) 
D(t) 
E(t) 



t rt 



JO 
t r t 



JO 
t rt 



o Jo 



(v s (0v s (0Uv c (0v c (0) c - (y s (0)s(v s (c)) s (K(0)e(K(C)>. 
(v s (0)s(v s (0) s [(v c (Ov c (0} e - (K(0)e(K(C))c]}« 
(v s (Ov s (0) s - (v s (0) s (v s (0) s ] (K(0)c(K(C))ci«- 



(l|e,s« 



|l)c 



(21) 



The above correlation integrals (|21|) are starting point for our theoretical investigations. 
In certain situations they can be simplified even further. We will study four different 
regimes in which simplification is possible: (i) mixing regime (corresponding to chaotic 
classical dynamics in both subspaces) in which arbitrary correlation functions appearing 
in fl2~T|) decay to and their integrals thus grow as oc t, (ii) regular regime in which due 
to absence of mixing the whole correlation integral C(t) (or D(t), or E(t)) grows as 
oc t 2 , and two regimes in which we have a separation of time scales, with the time scale 
of the environment being much shorter than the time scale of the central system. In this 
case we will work out two different regimes depending on the mixing property of the 
environment, namely (iii) "fast mixing" environment where environmental correlations 
decay; and (iv) "fast regular" environment where the environmental correlation function 
has a non- vanishing time average value. The decay of fidelity and purity-fidelity in 
mixing and regular regimes has already been extensively discussed || pi], |Tl|, |12|. In 
this work we do not only discuss a new quantity, namely the reduced fidelity, but also 
the unperturbed dynamics is now separable and thus not ergodic on the total space. 

It is interesting to notice that the general inequality (|13|) is clearly satisfied in our 
linear response results since the functions E(t) and D(t) are written in terms of sums of 
non- negative real numbers [eqs. ([18], 19)] and are hence themselves always non-negative. 

We want to stress that all the results of this section can be directly translated to 
the discrete time case of quantum maps (kicked quantum systems) by simply treating t 
as integer variable and replacing all the integrals by sums, J * — > X^ -1 - 



4. Numerical experiment: two coupled kicked tops 

We will now illuminate and demonstrate our theoretical predictions with a numerical 
example of two coupled kicked tops. In addition, numerical simulations will provide 
us with some insight about the asymptotic behaviour for long times and small 
fidelity/purity. Note that here the time is discrete (integer) and measured in the number 
of kicking periods (steps). A single kicked top |L6| has a unitary one step (Floquet) 
propagator 

U(a, 7) = exp (-hJ y ) exp (-ia J 2 J2 J). (22) 

The level of chaoticity (e.g. the rate of mixing) of a single top can be varied by varying 
a and the time scale can be influenced by changing the angle 7. In order to be able 
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to study the reduced fidelity and the purity we couple two kicked tops, where one top 
will act as an environment an the other as a central system. Two coupled kicked tops 



14, 10] have the following unitary one step propagator 

U s := U s (a s , %)U e (a e , j e ) exp (-iS -V/K), (23) 

with C/ Sje being propagators for a single kicked top of a system or environment. The 
unperturbed propagator U := Ug=o is simply obtained by putting 5 = into above 
expression (123|) . Perturbation V = V s Cg) V e will be chosen to be of two different forms: 
for regimes (i),(ii) and (iii) both V^ e = J z /J have the same form, whereas in the regime 
(iv) we will take V s = J z /J and V e = J 2 / J 2 . The reason for chosing a different form 
of perturbation in case (iv) is that we there want the environmential time correlation 
function to have a non- vanishing time-average in order to yield generic results. In other 
cases (i-iii) the precise form of the force and position operators is irrelevant, so V = J z /J 
provides the simplest choice. The Planck constant is determined by J as h = 1/J, so 
that the semiclassical limit implies J — > oo. The initial condition will always be a direct 
product of coherent states for both tops, |1, 1) = |y> SJ $s) (g> l^e,^), with the expansion 
of SU(2) coherent state |<£>, i?) into eigenbasis \m) of J z being 

J / O J \ 1/2 

1^,0} = Y ) cos(^/2) J+m sin(^/2) J - m e- fm ^|m). (24) 

Let us now work out the details of the fidelity, reduced fidelity and purity decay in all 
four mentioned regimes. 



4-1. Mixing regime 

If the combined correlation function decays to zero sufficiently fast, so that its time 
integral converges and C(t) oc t (|TTD, we can define a kind of transport coefficient 

o := lim C{t)/2t. (25) 

t — >oo 

Linear response formalism thus gives us initial linear decrease of fidelity. On the other 
hand, the coefficient D(t)+E(t), occurring in the expression for purity, is small compared 
to C(t) 0, namely {D(t) + E{t)}/C{t) ~ 1/N S + l/N c (N s>e are dimensions of H s>e ) 
and the same holds for each individual E(t) or D(t). Therefore, up to semiclassically 
vanishing correction, all the three quantities are expected to decay on a same time scale. 
What is more, if also multi-time-correlation functions fall off fast |], [10], [TTf] , then the 
shape of decay for longer times is exponential with the exponent given by the linear 
response formula, 

F\t) = Fl{t) = I{t) = exp (-2t/r m ), r m = JjL-. (26) 

This formula is expected to be valid for times longer than the Ehrenfest time or classical 
mixing time, depending on the initial coherent or random state, respectively, and for 
sufficiently small S such that the value of fidelity/purity is still far above the saturation 
for the above mentioned time scale (see |I| for a detailed discussion). Therefore, if the 
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central system and the environment are both separately mixing (by mixing we mean 
decay of all time correlation functions) the fidelity, reduced fidelity and purity decay in 
the same way. This means that the decay of purity, decay of reduced fidelity and decay 
of fidelity have the same principal mechanism, that is "corruption" of the state vector 
as a whole which gives dominant contribution over more subtle issues of destruction 
of coherence. Note that stronger mixing (usually connected to stronger chaoticity of a 
classical system) means smaller correlation integral a and thus slower decay of purity, 
fidelity and reduced fidelity. This is totally opposite to the very short time behaviour 
described by Refe.g, @, Q. 



0.1 




o 



1e-05 



10 



100 



1000 



t 



Figure 1. Correlation sums (|2l]) in the mixing (chaotic) regime (see text) divided 
by 2i: C(t) (solid), C(t) - D(t) (long dashed), C(t) - D(t) - E(t) (short dashed). 
Horizontal (chain) line shows best fitting value of a j25|). 



For numerical verification of this result we chose the perturbation V^ e = J z / J and 
parameters a Sje = 30, 7 Sje = vr/2.1, and J = 1/h = 200 giving the effective size of the 
Planck constant. Coherent initial state is placed at i? Sje — tt/ tp Sje = ir/ \f2. In figure 
|l] we plot time dependence of correlation integrals occurring in expressions for fidelity 

]18|) , and purity (19), showing that the terms D{t) and E(t) are 



(llq), reduced fidelity 

really negligible as all the three quantities shown are practically equal, for times longer 
than the Ehrenfest time. Next we show in figure |] the decay of fidelity F(t), reduced 
fidelity F R (t) and purity I(t) for the same parameters and for the perturbation strength 
5 = 8- 10~ 4 . Clean exponential decay is observed in all three cases, on a time scale r m 
(p6|) given exactly by the lowest order linear response expression (|16D in terms of a (p5|) 
which is obtained from data of figure [l]. Exponential decay, of course, persists only up 



to the saturation value determined by a finite Hilbert space size [10 . 
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Figure 2. Decay of F 2 (t),F^(t) and I(t) (dotted curves) in the mixing (chaotic) 
regime. The solid line gives exponential decay ( |26| ) with r m calculated from a in figure 
[p. Horizontal chain lines give the saturation values of purity and reduced fidelity, 1 /200 
and 1/400 2 , respectively. 



4-2. Regular regime 

If the system is regular as a whole, then the time correlation function will generally not 
decay to zero, but will have a non-vanishing average value 

c F := lim C(t)/t 2 . (27) 

t — >oo 

Similar coefficients can be defined for the average of C(t) — D(t) x cr^ 2 occurring in the 
expansion of the reduced fidelity, and C(t) — D(t) —E(t) x cit 2 for the purity. Note that 
C(t) is proportional to H for coherent initial states. As shown before [TT|, the expression 



D(t) + E(t) is almost equal to C(t) for coherent initial states in the regime of regular 
- integrable classical dynamics. The expression C(t) — D(t) — E(t) occurring in the 
formula for purity (|l~9|) is therefore of the order h 2 and the decay time scale for purity 
decay is h independent. This cancellation in the leading order in h happens due to both 
E(t) and D(t) terms and the reduced fidelity will therefore decay on approximately the 
same time scale as the fidelity. For coherent initial states one can show [|K], |Tl| that the 
shape of decay is a Gaussian 

F(t)=exp(-(t/r r ) 2 ), r r = -^=, (28) 

and with similar expression, only with cr replacing cp, for reduced fidelity. Remember 
that for coherent initial states cf,r oc h and therefore r r oc yh. The time scale for decay 
of purity is again given by an analogous expression, namely h/8\/c\, with cj oc h 2 , but 
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as we will see in the numerical simulations and discuss later, its long time behaviour is 
not a Gaussian but has an algebraic tail instead. 

For the numerical demonstration we take two kicked tops with V^ e = J z / J, J = 200, 
7s,e — tt/2.1 and a Sje = in order to have regular dynamics in both subspaces. Initial 
coherent state is placed at ("&,(p) B = (k/VS, n/V2) and (i?, (p) e = (7r/\/3, 37r/v7). In 




1e-06 




1000 



Figure 3. Correlation sums ( pl|) in the regular regime (solid curves). Chain lines 
show the theoretical values of cf,r,i (|29|). See text for details. 



fact, the dynamics for 7 es = 7r/2 (approximating a 'more generic' value 7 = vr/2.1, 
and a = 0) is very simple, namely it is a pure rotation by an angle 7r/2, and 
all three correlation integrals (sums) can be calculated explicitly. Coherent state 
expectation values of the time averaged correlation functions are then easily evaluated 
using formulas for the expectation values of the lowest powers of J x , J y , J z . If we denote 
with n se = (x,y,z) B)e = (sin?? cosy?, sin sin 93, cos$) S)e the unit vectors in the direction 
of initial coherent states of the central system and the environment, then the results for 
c (for this special case of 7r/2 rotation with a = 0) are 



ci 



8J 
1 

8J 
1 



y s 2 -y c 2 -2(n s 



n e - y s y c 



+ 



n s ■ n n 



VsVc 



16J 2 



{Us ~ Vcf + (n s • n e 



+ 



VsVe 



16J 2 
,2' 



16 J 2 



(y s - Ve) 2 + (n s ■ n c -y s y c 



y s y e f 



(29) 



Note that all c's are expressed in terms of the invariants of motions as they should be 
by definition (however, this does not mean that the perturbation itself is an invariant of 
motion). Here one can also explicitly see the cancellation of terms for purity, namely the 
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cp and cr are proportional to H, while c\ is proportional to h 2 . We should note that for 
special positions of initial packets, the average correlator may vanish c = 0, and there 
the decay may be much slower and not Gaussian at all. The zeroes of c's are therefore 
very special points denoting wave packets that are very stable against perturbations. 
For fidelity and reduced fidelity this slow decay can give rise to a power law decay of an 
average fidelity (fidelity averaged over the whole phase space). 

In figure |] we can see that the three theoretical coefficients c's ( p9|) for 7 s e = n/2 
approximate very well the numerical calculation of correlation sums in the case of 
7s,e = vr/2.1. In figure |] we show the decay of fidelity, reduced fidelity and purity 




1 10 100 1000 10000 

t 

Figure 4. Decay of F 2 (t), F^(t) and I(t) (solid curves) in a regular regime (see text). 
Dashed curve is a Gaussian with the exponent given by c\. Note that theoretical 
Gaussians for F(t) and Fn(t) are indistinguishable from the data! Dotted line has a 
slope -2, and the horizontal chain line gives the saturation value of I{t) at « 1/130. 



for 5 = 5 ■ 10 . The decay of fidelity and reduced fidelity is Gaussian fl28|) with the 
decay times r r given very accurately by the theoretically calculated cf and cr shown 
with solid curves in figure [| The decay of purity on the other hand is not Gaussian for 
long times. Of course, it decays quadratically as given by a linear response formula ([191) 
for times short enough that the purity is close to 1, but for larger times it decays with a 
power-law. Still, up to a constant numerical factor (independent of S, % etc) , the effective 
purity decay time is given by a linear response formula for r r (^). This can be observed 
in figure |5] where we show numerically calculated values of r r based on 1/e m 0.37 level 
fidelity (purity) together with the theoretical prediction on a basis of average correlation 
function (c). Since the purity decay does not follow a Gaussian model the numerical and 
theoretical scale deviate by a numerical factor, which is however constant (independent 
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Figure 5. Times r at which F 2 (t),F^ (t), I(t) fall to level 0.37 for different 6 in a 
regular regime. Symbols: pluses, circles, and triangles, are numerics whereas lines: 
solid, dashed, dotted, give theoretical dependence of r r , for F 2 (t),F^(t), and I(t), 
respectively. All parameters are the same as for fig. and |L except here: J = 100. 



of h and 5). For larger perturbation strength S, the decay rates become less and less 
sensitive to perturbations in agreement with the expected saturation 0. We stress 
that for short times, namely for 5t -C 1, the purity of a Gaussian packet for regular 
dynamics can be written as an algebraic function in terms of determinants [|ll|]. This 
dependence is completely independent of % and is a function of the scaling variable 
5t alone. Furthermore, we expect periodic revivals of purity with a classical beating 
frequency v ~ 1/5, again independent of h. Based on numerical results (extending the 
analytical results for longer times) we conjecture that the overall dependence of I(t) 
is h (or J) independent and its functional form can be described in terms of a scaling 
variable St, namely I = f(St). This is demonstrated in figure |6]. 

The regular regime discussed above has practical importance for the emergence of 
the macro world [T7j . If we have a macroscopic superposition of states the decoherence 
time is going to be smaller than any dynamical time scale in a system and one is trivially 
in a regular regime with non-decaying correlations. 



5. Separation of time scales: Fast environment 

The expressions for C(t), D(t) and E(t) can be further simplified if the decay time scale 
of the environmental correlations (V e (t)V e (t')) e is much smaller than the time scale of 
the systems' correlations (V s (t)V s (t')) s . The time averaging over the fast environmental 
part of the perturbation V e can be performed in this case. Regarding the environmental 
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10 20 30 40 50 60 70 

5t 

Figure 6. Decay of purity I(t) in a regular regime for J — 200 (full curve), J = 100 
(dashed) and J = 50 (dotted) and the same parameters as in fig. |]. Note that I(t) 
is almost fully independent of h, and that curves for other values of S are completelly 
overlaping with the existing curves. 

correlation function two extreme situations are possible. Namely, the correlations of 
the environment decay ( "fast mixing environment" ) so that we have a finite integral of 
environmental correlation function, or the correlations of the environment do not decay 
( "fast regular environment" ) and we have generically nonvanishing average correlation 
function of the environment. 

5.1. Fast mixing environment 

The situation, when the time scale t e on which correlation functions for the environment 
decay is much smaller than the time scale t s of the central system, is of considerable 
physical interest. This includes various "brownian" like baths, where correlation times 
are smaller than the dynamical times of the system in question. The expressions for 
C(t), D{t) and E{t) (|2l|) can be significantly simplified in such a situation. We will 
furthermore assume (V c ) c = 0, with (A) = lim^oo t^ 1 J * (A(£))dt; denoting the time 
average, which is true if we are in an equilibrium situation (average "force" V e vanishes). 
The integration over the fast variable V e in (|2l|) can then be carried out and we get 

Cit) = 2ta c W) s 

C(t)-D(t) =2ta e {(V s \-(V s )\} 

C (t) - D(t) - E(t) = 2ta c ~ (f s j , (30) 
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with 

a e := lim (S c 2 (t)) c /2t, and E e (t) = f K(£R, (31) 

being the integral of autocorrelation function for the environmental part of perturbation 
V e alone. From expressions (|3"0|) and the lowest order expansions ( |I1^|T8| , |I§D we can 
see that the decay time scale depends only on the time average diagonal correlations 
of central system (V s (t) 2 ) and not on the full correlation function. This is a simple 
consequence of the separation of time scales and means that the decay of all the three 
stability measures does not depend on the dynamics of the central system (e.g. being 
mixing (chaotic) or regular.) Furthermore, the reduced fidelity F&(t) and the purity 
I(t) will decay on the same time scale (|3~0l), meaning that the decay of reduced fidelity is 
predominantly caused by the loss of coherence, i.e. entanglement between the two factor 
spaces. This means that the reduced fidelity, which is a property of echo dynamics of 
comparison of two slightly different hamiltonian evolutions, is equivalent to the decay 
of purity or growth of linear entropy of an individual weakly coupled system. 

If the initial state of a central system p s (0) is a Gaussian wave packet (coherent 
state) then the dispersion of (V B ) S — (K) 2 S is by a factor of order 1/h smaller than (V^) s . 
Thus for coherent initial states of a central system, no matter what the initial state of the 
environment, the Fn(t) and I(t) are going to decay on a 1/h times longer time scale than 
F(t). We have therefore reached a general conclusion based on very weak assumptions 
of chaotic fast environment, namely that the coherent states are most robust against 
decoherence (provided t c <C £ s ), and that decoherence takes place in times longer than 
the correlation time of environment t e <C tdec- If decoherence is even faster than the time 
scale of the environment, as is the case for macroscopic superpositions, then formulas 
( ftOf) are not valid any more as one is effectivelly in a regular regime of the previous 
section. Decoherence time is then independent not just of systems dynamics but also of 



environmental dynamics characterized by o c (see ||17||). 

In the regime of fast chaotic environment one can immediately derive a master 
equation for a reduced density matrix of a central system Jl8|, |19|| . We take partial trace 
over the environment of expansion for pMs{t) (|5|, p!5|) and write it for a small time step At. 
This time step At must be larger than the correlation time t e of the environment and at 
the same time smaller than the correlation time t s of the system. For the environmental 
part of the correlation function we assume fast exponential decay (particular exponential 
form is not essential) which is independent of the state p 

tr e {V e (t)V e (t')p} — > ^exp{-\t-t'\/t c }ti cP . (32) 

te 

Assuming the perturbation to be a product V(t) = V a (t) <8> V e (t) and the average "force" 
tr c (V c (t)p) to vanish together with the exponential decay of environmental correlations 
of the form (|32|) for arbitrary state, yields a master equation for the reduced density 
matrix p Ms (t) := tr e p MS (t) 

PMs(t) = --£\y B {t), [V B (t),p Us (t)}}- (33) 
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Figure 7. Various correlation sums from formulas (|3C|) in a fast chaotic regime (solid 
curves, as indicated in the figure). Chain lines indicate corresponding theoretical time 
averages. For details see text. 
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Figure 8. Decay of F 2 (t),F^(t) and I(t) for fast chaotic environment. Dashed line is 
exponential with the exponent given by the values of <j e and (V r s 2 ) s ( pcf ) and two dotted 
lines have slopes —2 and — 1. 



For numerical demonstration we choose: V Sie — Jz/J> J — 200, 5 = 1.5 • 10 3 , 
coherent initial state at <^) s>e = (n/y/3, tt/V2) and parameters a s = 0, 7 S = 7r/50 for 
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Figure 9. Times r at which F 2 (t), F£ (t),I(t) fall to level 0.37 for different 6 and fast 
chaotic environment. Symbols give numerics and lines give theoretical dependence of 
r (same as in fig. ||). All is for J = 100. 



the central system and a e = 30, 7 e = tt/2.1 for the environment. Actually, we could 
take any value of a s and would get qualitatively similar results. The only advantage of 
using regular central dynamics a s = is that it is then possible to explicitly calculate 
averages (K?) s and (V^. 2 ) s — (V S ) 2 S . Namely, if a s = and 7 S C 1 we get 

(W)s-(^Fs = ^(l + 2/s)- (34) 

The values of these two quantities for our initial condition are shown in figure [5] with 
two dotted lines (by pure coincidence we have a e ~ (V 2 (t)) s ), together with numerically 
calculated time dependent (not yet averaged) (V 2 (t)) s and {V 2 (t)) s — {V s (t)) 2 for our 
choice of 7 S = tt/50. This time dependent values oscillate on a time scale ~ 50, which 
is much longer than the time ~ 10 in which cr c (|3~TD converges and so the assumption 
t e <C t s is justifiable. The value of all these three quantities is then used in linear 
response formulas (^0|) to give us time scales on which F, Fr and I decay. The results 
are shown in figure |8]. We can see that the fidelity again decays exponentially, but the 
reduced fidelity and purity have a power-law like tails. They decay on a time scale still 
roughly given by the lowest order expansions (|30]) and the values of o~ c (numerical from 
fig. 0) and (V^, 2 ) s — (V S ) 2 S (theoretical expression (0)) which can be seen in figure []. The 
same general conclussion again holds: the more chaotic the environment is (smaller o~ c ), 
the slower the decay of all three quantities. Purity and reduced fidelity both decay on a 
1/h longer time scale than the fidelity in accordance with expressions (0) for coherent 
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initial states. 

5.2. Fast regular environment 

Here we will explore perhaps a less physical situation of a regular environmental 
dynamics. For a regular environment the double integral of environment correlations 
grows as oc t 2 due to non-decay (plateau) of correlation function and we can define the 
average correlation function 

c e := lim (S c 2 (t)) e /t 2 . (35) 

t — >oo 

If in addition the correlations of the system also do not decay then the correlation sum 
of the total system will grow as oc t 2 which is a regular regime discussed previously. 
Here, we will focus on a different situation where the integral of system's correlation 
function converges, that is C s (t) oc t, i.e. the central dynamics is mixing (chaotic). We 
will additionally assume the average "position" V s to be zero (V a ) = 0. The transport 
coefficient of a system a s is then 

cr s := lim (S s 2 (t)) s /2t, E a (t) = f V*(£)d£. (36) 

The expressions for C(t), D(t) and E(t) (|2T|) for the present case can be simplified to 
C(t) = 2ta s c c 

C(t) - D{t) = 2ta s c e 

C(t) - D{t) - E(t) = 2ta s {c c - ml} • (37) 

An important thing we notice immediately is that now the reduced fidelity F&(t) decays 
on the same time scale as fidelity F(t). This must be contrasted to the case of a fast 
mixing environment (P0|), where Fn(t) decayed on the same time scale as purity. If the 
initial state of the environment p e (0) is a coherent state, then the purity will decay on 
a 1/h times longer time scale as fidelity and reduced fidelity. On the other hand, for a 
random initial state of the environment, the average force vanishes (V e ) e = 0, and then 
all the three quantities decay on the same time scale. 

For the purpose of numerical experiment we chose now V s = J z /J, and V e = J 2 / J 2 
in order to have a less trivial situation of non-vanishing average force. The initial 
condition is again (f) s ,e = (k/V3,it/\/2) and parameters are J = 200, ct s = 30, 
7 S = 7i77 and cx e = 0, 7 e = it 12.1 and perturbation strength 5 = 6- 10~ 4 . By choosing 

2 

explicitly solvable case a e = we can calculate c e and c c — (V e ) e , say for the simple case 
of n / 2 rotation, 7 e = tt/2, where we obtain 

ce = \(l ~ Vlf + ^(-3y c 4 + 2y 2 + 1) + 0(1/ J 2 ) 

ce - We = - Vt) + J^p( n yt - n Ve + 2 ) + (VA (38) 

The values of this coefficients are shown in figure |TD| (lower two dotted lines). They 
agree nicely with numerics for 7 e = 7r/2.1. In the figure |Tl| we can observe exponential 
decay of fidelity and reduced fidelity on the same time scale (two curves almost overlap) 



Fidelity and Purity Decay in Weakly Coupled Composite Systems 



18 



10 



0.1 



0.01 



0.001 



0.0001 



I I I 1 1 1 




/<E s 2 >/2t 










V 2 ft 2 

<E e >/t 










C e (t)/t 2 : 


1 1 1 . . 




10 


100 


1000 



t 

Figure 10. Correlation sums occurring in (^) (solid curves). Top chain line gives best 
fit for cr s and two lower chain lines give theoretical time averaged correlation functions 
for the environment (Bq). All is for a fast regular regime. See text for details. 
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Figure 11. Decay of F 2 (t), F^(t) and I(t) for fast regular environment. Dashed line 
is an exponential with the exponent given by a product of cr 8 and c c (|37]). Straight 
dotted line has a slope —1. See text for details. 



and decay of purity on a 1/h longer time scale. For longer times the purity decay is 
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Figure 12. Times r at which F 2 (t) (pluses), F#(t) (circles), I(t) (triangles) fall to 
level 0.37, and the times r when I(t) falls to 0.99 (squares) for varying S and with fast 
regular environment. Symbols are numerics and lines give theoretical dependences of 
r. Everything is for J = 100. 



again algebraic. In figure 12 we show dependence of decay times on 5. The dependence 



for purity is quite interesting. If one looks at the time the purity falls to 0.99 one has 
the agreement with linear response (by definition). But if one looks at the purity level 
0.37, they don't agree as well, meaning that the nature (shape) of purity decay may 
change (not only the scale) as one vary 5 or %. On the other hand, this may also be 
simply a finite size effect due to finite Hilbert space dimensions. 

6. Discussion and conclusions 

In this paper we have analyzed the stability of unitary time evolution of composite 
systems under the weak coupling between subsystems. This is a natural (unitary) model 
for dissipation and decoherence in quantum mechanics where one subsystem plays the 
role of the central system and the other plays the role of the environment. But this is 
not the only possible application of the above ideas. One may also be interested in the 
dynamical effects of weak coupling between two controllable parts of the system, like 
e.g. the atom and the electromagnetic cavity |2(| . 



We have analyzed three different quantities which are treated on a similar 
theoretical footing but which have different physical interpretations. The first two, 
namely the fidelity and the reduced fidelity, refer to the case of echo dynamics where 
the forward evolution is generated by the uncoupled system while the backward evolution 
is generated by the weakly coupled system. The third one, namely the purity, refers 
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to the growth of linear entropy or growth of entanglement between the two subsystems 
during the course of weakly coupled forward evolution. First we have shown a rigorous 
inequality between the three quantities which may be useful to provide various bounds. 
Then we have developed a linear response theory which predicts time scales for the 
three quantities in terms of time correlation functions of the perturbation in each of the 
subsystems. Thus we have been able to classify all different behaviours with respect 
to regularity or chaoticity (as defined here by the mixing property) of each of the 
subsystems. The general conclusion is again, consistent with |22|| , that strong chaos 
stabilizes quantum dynamics with respect to the intersystem coupling, and that strong 
chaos decreases the rate of entanglement (or linear entropy) growth. In particular, if 
the characteristic time scales of correlation decay in two subsystems are well separated, 
then we can integrate over the correlation functions of the fast subsystem and obtain 
expressions which are independent of the nature of dynamics in the slow subsystem. 
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